Phase transition in two dimensions far from thermal equilibrium 



Patrick Dillmann, 1 Georg Maret, 1 and Peter Keim 1 

1 University of Konstanz, D-78457 Konstanz, Germany 
(Dated: March 28, 2013) 

A two-dimensional monolayer of monodisperse colloids is quenched rapidly below the melting 
temperature into the crystalline phase and the nature of the transition is investigated. Analyzing 
the local order as function of time we find shortly after the quench a bimodal distribution of highly 
ordered (crystalline) and fluidlike particles, indicating that the KTHNY-scenario known from 2D 
equilibrium transitions does not hold. Critical nucleation theory can not describe the solidification 
scenario found in our experiment, either, as we do not find a lag time nor a critical nucleus size. 
The average size of nuclei grows monotonically as expected, but surprisingly, comparing the size of 
an individual 'nucleus' in consecutive timesteps, shrinking is always more probable than growing. 
This indicates a strongly asymmetric nucleation behaviour. Nevertheless, the averaged dynamics of 
the degree of crystallinity is found to be in very good agreement with the heuristic Johnson-Mehl- 
Avrami-Kolmogorov(JMAK)-scenario. 
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Crystal nucleation is one of the basic physical phe- 
nomena that govern the growth and the structure of any 
crystalline state in three dimensions. Therefore this topic 
[THIS] is one of the most important and interesting areas 
of condensed matter physics. Although it is rather easy 
to control the crystal nucleation process empirically, un- 
derstanding this phenomenon at the fundamental level 
has been a great challenge for many decades. The most 
widely used concept, the Classical Nucleation Theory 
(CNT), attempts to describe this complex process by 
treating the nuclei as compact spheres that grow if their 
size exceeds a critical value [THlj. This critical value is 
determined by a maximum in the free energy, where the 
energy gain in the volume ~ r 3 overcompensates the sur- 
face energy cost ~ r 2 of the nucleus. Colloidal systems, 
where micrometer sized particles are dispersed in a sol- 
vent, allow to investigate such nucleation on a single par- 
ticle level. Depending on interaction strength and volume 
fraction of the colloids in the solvent, different thermo- 
dynamic phases like crystalline, fluid, gaslike and (in two 
dimensions) hexatic phases can be identified. Colloids 
are small enough to undergo Brownian motion and can 
be described as a statistical ensemble but they are large 
enough to be monitored with optical methods like light 
scattering or (confocal) video microscopy. But not only 
the structure can be monitored on an 'atomistic' scale. 
Since the particle dynamic scales inversely with the size of 
the particles, structural rearrangements like phase transi- 
tions can be monitored on single particle level on almost 
all relevant time scales. 

However, experiments on hard sphere colloidal systems 
[HI E] , where individual particle interactions are suffi- 
ciently simple to be modeled theoretically differ from the 
predictions of critical nucleation theory. In more complex 
systems, the CNT calculation underestimates the nucle- 
ation rate by many orders of magnitude [16] . It is obvious 
that the simple picture of a spherical seed has to be ex- 



tended. In computer simulation studies of hard sphere 
colloidal systems in 3D [17J . and for phase-field-crystal 
calculations in 2D [18], a precursor-mediated two-step 
crystallization was observed. Dense amorphous clusters 
are formed as the first step and then those clusters act 
like precursors for the nucleation of well-ordered crys- 
tallites. In experiments in 2D with effective attractive 
interactions between particles such precursors are found 
[TO1 |2T)] . too, but might be explained with sublimation 
rather than nucleation. Dusty plasma experiments in 2D 
showed a power-law increase of average crystallite do- 
main size as measured by bond-order correlation length 
[21) and microheterogeneity in structural evolution [22j . 

It is well known that dimensionality affects phase tran- 
sitions. Whereas nucleation of supercooled fluids in three 
dimension is a first order phase transition including phase 
equilibrium, the melting process in a two-dimensional 
system in equilibrium has been a matter of debate for 
long. While Kosterlitz and Thouless proposed a contin- 
uous transition induced by topological defects like dis- 
locations at the melting temperature T m [331 123 > Chui 
later argued for a first order transition similar to the case 
in three dimensional systems |25| . Here, the symmetry 
is restored in the high temperature phase by the spon- 
taneous generation of grain boundaries. The presence 
of a hexatic phase, an intermediate phase between the 
isotropic liquid and the crystalline state, was predicted 
by Nelson and Halperin [25J [57] , as extension of the the- 
ory by Kosterlitz and Thouless. The hexatic phase is 
an anisotropic fluid and a second type of topological de- 
fects, namely disclinations, are needed to restore isotropy 
at a higher temperature Tj. Young independently cal- 
culated the critical exponents at the melting tempera- 
ture T m [2H]. As colloidal systems are excellent model 
systems to observe phase transitions, experiments have 
been done on 2D or quasi-2D colloidal systems using op- 
tical microscopy with video imaging 29 33J to verify the 
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predicted two step transition and the existence of a hex- 
atic phase. In [3H [35] we were able to investigate the 
microscopic mechanism of the so-called KTHNY-theory 
based on defect unbinding. The critical exponents of di- 
verging length scales are in agreement with renormaliza- 
tion group calculations. In heating and cooling cycles we 
found the KTHNY-picture to hold as melting and freez- 
ing scenario without hysteresis of the transition points, 
but care has to be taken that the system is always in 
thermal equilibrium. We did not observe a phase equi- 
librium or other first order signatures and accordingly we 
identified the transitions to be continuous. 

Here, we present a study of solidification far from ther- 
mal equilibrium. If the system is quenched rapidly at 
ultra-short timescales from the isotropic fluid phase into 
the crystalline phase the scenario alters completely: we 
did NOT observe the two distinct phase transitions and 
the intermediate hexatic phase [36] . If quenched rapidly, 
the system solidifies into a poly-crystalline sample. We 
investigate this phase transition far from thermal equilib- 
rium on single particle level and show that critical nucle- 
ation theory does not describe the phenomena observed. 

The 2D colloidal monolayer is realized out as follows. A 
droplet of a suspension composed of polystyrene spheres 
dispersed in water is suspended by surface tension in a 
top sealed cylindrical hole (0 = 6 mm) of a glass plate. 
The particles are 4.5 fim in diameter and have a mass 
density of 1.5 g/cm 3 leading to sedimentation. The large 
density is due to the fact that the polystyrene beads are 
doped with iron oxide nanoparticles leading to superpara- 
magnetic behavior. After sedimentation particles are ar- 
ranged in a monolayer at the planar and horizontal water- 
air interface of the droplet and hence form an ideal 2D 
system. Particles are small enough to undergo Brownian 
motion but large enough to be monitored with video- 
microscopy. An external magnetic field H perpendicular 
to the water-air interface induces a magnetic moment in 
each bead (parallel to the applied field) leading to repul- 
sive dipole-dipole interaction between all particles. We 
use the dimensionless control parameter F to character- 
ize this interaction strength; it is given by the ratio of 
dipolar magnetic to thermal energy 
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and thus can be regarded as an inverse system temper- 
ature. The state of the system in thermal equilibrium - 
liquid, hexatic or solid - is solely defined by the strength 
of the magnetic field H since the temperature T, the 
2D particle density p and the magnetic susceptibility 
per bead \ are kept constant experimentally. In these 
units the melting temperature (crystalline - hexatic) is at 
about r = 60 and the transition from hexatic to isotropic 
at about T = 57. Based on an equilibrated liquid system 
(r « 15) [37j we initiate a temperature jump with cooling 
rates up to dT/dt ~ 10 4 s _1 into the crystalline region of 
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FIG. 1. The color code measures the frequency of locally 
ordered particles in the m-n-plane 220 s after a (inverse) tem- 
perature quench from F = 13 (deep in the fluid phase) to 
r = 63 in the crystalline stage. A bimodal distribution of 
crystalline-like particles with high local order (upper right 
corner) and fluid-like particles with low local order (lower left 
corner) is found. 



the phase diagram Tm > 60 [3H [35] . This temperature 
quench triggers a crystallization process in the system 
where the time scale of cooling is 10 5 faster compared 
to the intrinsic Brownian timescale Tg = 50 sec, i.e. the 
time particles need to diffuse the distance of their own di- 
ameter. An elaborate description of the experiment can 
be found in [38]. The field of view is 1160 x 865 pm 2 in 
size and contains about 9000 colloids. Each temperature 
quench to the same final value of the control parameter 
IV is repeated at least ten times with sufficient equili- 
bration times in between and all analyzed quantities are 
averaged by the number of repetitions of the quenches. 

To analyze the degree of local order after a temperature 
quench we use a criterion based on the magnitude me k = 



\ipk\ °f the local bond order parameter 
local bond order field is given by 
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Here the sum runs over the Nj next neighbors of the 
particle k at position r and 9jk is the angle between a 
fixed reference axis and the bond of the particle k and 
its neighbor j. For particle k, m§ k is close to one if the 
neighboring particles are sitting on a hexagon and close 
to zero if the neighboring particles are not sixfold. One 
can sharpen this measure of local sixfold symmetry by 
taking the magnitude of the projection of ipk 
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to the mean local orientation field given by the nearest 
neighbours. It takes the second nearest neighbors into 
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FIG. 2. Average size of the nuclei (in units of the voronoi- 
cell of individual particles) for five different quench depth as 
function of time. The nonmonotonic behavior of the deepest 
quench is an artefact, since (large) crystallites which exceeds 
the field of view are ignored. Note, that nucleation starts 
immediately after the quench. 



account and determines how the orientation of particle k 
fits into the orientation of its neighboring particles. Fig. 
[T] shows the probability distribution in the mQ-n§-p\&iie. 
The color code measures how often a particle with high 
or low local order is found. Compared to the melting 
scenario in thermal equilibrium, where only a unimodal 
distribution is found [33 [40] , one can clearly identify a 
bimodal distribution compatible with a nucleation sce- 
nario. To investigate the time dependence of growing 
grains, we define a particle to be part of a crystal clus- 
ter if the following three conditions are fulfilled for the 
particle itself and at least one nearest neighbour: a) The 
magnitude of the local bond order field mg fc must exceed 
0.6 for both neighboring particles, b) The bond length 
deviation A\lki\ of neighboring particles k and I is less 
than 10% of the average particle distance l a . c) The 
variation in bond orientation AO — \ipi- — of neigh- 
boring particles k and I must be less than 2.3° in real 
space (less than 14° in sixfolded space). The threshold 
values are denned in comparison with mono-crystalline 
and isotropic fluid phases under equilibrium conditions 
and the qualitative results are robust to (moderate) vari- 
ations (about ±20%) of the threshold values. Simply 
connected domains of particles which fulfill all three cri- 
teria are merged to a crystalline cluster. Comparing our 
results with [17H20j . no amorphous precursors with in- 
creased density are found in our experiment. Analyzing 
crystalline and fluidlike domains the 2D density measured 
by Voronoi-cells is identical within the error bar. Figure 
[2] shows the average size of the crystalline cluster up to 
10 min after the quench. As expected, the average size of 
the nuclei grows monotonically as function of time, for all 
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FIG. 3. The shrinking probability after a quench (here 
Tf = 93.5 is shown) for a crystalline cluster composed of M 
particles is always larger than the growing probability. Unlike 
in 3D systems one can not determine a critical nucleus size on 
single particle level and hence 2D systems out of equilibrium 
can not be described with critical nucleation theory. Solidifi- 
cation is strongly affected by order fluctuations and shrinking 
and growing appear to be quite asymmetric. 



quench depths ranging from quenches close to the melt- 
ing temperature T m to quenches deep in the crystalline 
phase. We checked that the maximum in size for the 
quench to T = 141 is an artefact since only crystallites 
being completely in the field of view are counted. For 
this deepest quench some crystallites grow beyond the 
field of view after about 5 min and are ignored when de- 
termining the mean size. It is important to note that the 
nucleation starts immediately after the quench such that 
no lag time of critical nucleation is detectable. Addition- 
ally, with the parameters given above about 10% of the 
particles are identified as crystalline before the quench, 
deep in the fluid phase at T = 13, 7. Such a precritical 
nucleus consists in average of 2.7 particles and the life- 
time is a few seconds. The large amount (compared to 
3D systems) of precritical nuclei is due to the fact that 
the preferred local symmetry in 2D is sixfold, both in the 
crystalline as well as fluid phase. (If the criteria for crys- 
talline particles are tightened, the amount of precritical 
nuclei decreases but then a monocrystal in equilibrium is 
not detected as such any more). This already indicates 
that the formation of a nucleus in 2D after a tempera- 
ture quench is not a rare event and the critical nucleation 
barrier is very low. In other words, the size of a critical 
nucleus consists of two particles, being the smallest pos- 
sible cluster. In this case one would expect all individual 
nuclei > 2 particles to grow after the quench which is sur- 
prisingly not the case. In Figure [3] the difference of grow- 
ing p g and shrinking probability p s of nuclei is plotted 
for different time windows after a quench to Tp = 93.5 
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where M is the size of the nuclei given by the number 
of the particles which belong to the cluster. The proba- 
bility for shrinking and growing for individually labeled 
clusters of size M is given by p g / s — N 9 / s /Nm where Nm 
is the number of clusters of size M and N g or N s are the 
numbers of those clusters which are larger M T > M or 
smaller M T < M after a finite waiting time r. 

For all time windows and nucleus sizes the probabil- 
ity of shrinking is larger than the probability of growing. 
This is true for all quench depths performed, close to 
melting (Tp — 61) as well as deep into the crystalline 
phase (Tp = 141, not shown here). The waiting time 
was varied over almost two decades from 0.1 — 10 sec and 
no qualitative differences were found. Obviously shrink- 
ing (often, but small areas) and growing (less often, but 
larger areas) is quite asymmetric and the nucleation pro- 
cess far from equilibrium is strongly dominated by fluc- 
tuations. Practically all grains which are present during 
the quench disappear again. Two mechanisms appear 
to be responsible for this. First, small grains dissolve 
completely into fluidlike particles and second, grains of 
intermediate size fluctuate strongly in orientation such 
that they disrupt into smaller subdomains. Such subdo- 
mains frequently merge again but then the net contri- 
bution to the growing probability is typically zero, since 
one grain has grown but the other (s) has/have 'disap- 
peared'. Interestingly the lifetime of crystalline clusters 
after the quench is again a few seconds and comparable 
to the lifetime of 'subcritical nuclei' before the quench. 

To further analyze the crystallization dynamics we do 
not measure the size of individual grains but determine 
the overall fraction of crystalline particles as function of 
time X(t) = NX N ^ W where N x ~ tal (t) is the number of 
crystalline particles and N(t) is the number of all parti- 
cles in the field of view at time t. Following a heuristic 
description of nucleation introduced by Johnson-Mehl- 
Avrami-Kolmogorov (JMAK-theory) the rate of crys- 
tallinity can be described with 
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if the following criteria are fulfilled: i) Nucleation occurs 
randomly and homogeneously over the entire untrans- 
formed portion of the material, ii) Growth occurs at the 
same rate in all directions, iii) The growth rate does not 
depend on the extent of transformation [41 - 43 . Whereas 
the first two criteria are reasonable, the last is complied 
since our system is not transport limited (the area den- 
sity of fluidlike and crystalline particles is the same) . 

Figure [4] shows the degree of crystalline particles as 
function of time. Solids lines are a fit to Eqn. [IJwhere the 
parameter to accounts for the 10% particles being already 
crystal-like in structure before the quench. Very good 
agreement is found and by taking twice the logarithm of 
Eqn. [| one can extract the Avrami-exponent (see inset). 
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FIG. 4. Fraction of crystalline particles for different quench 
depths. The solid lines are fits to Eqn. [4] The inset shows 
an Avrami-plot for Fe ~ 61.2 to determine the exponent n 
in Eqn. [4] from the slope of the curve. 



A: 



to [a] 



n 



61.2 
63.3 
69.8 
93.5 
141 



0.1 ± 10~ 7 
0.1 ± 10~ 5 
0.09 ± 10~ 7 
0.09 ± 10~ 6 
0.1 ± 10~ 4 



0.29 ± 10~ 7 
0.31 ±10~ 5 
0.35 ± 10~ 7 
0.4 ± 10~ 6 
0.45 ± 10~ 6 



0.09 ± 10~ 3 
0.13 ± 10~ 3 
0.28 ± 10~ 3 
1.27 ± 10" 3 
1.7 ± 10~ 4 

Whereas the parameter k is independent of the quench 
depth, the Avrami-exponent increases monotonically 
with the quench depth but is always smaller than unity. 

In conclusion, we found that neither the KTHNY- 
scenario nor critical nucleation does describe the 
phenomena observed for this 2D non-equilibrium 
scenario [44j . There is evidently no nucleation bar- 
rier since a) nucleations starts immediately after the 
quench and no lag time is observable, b) the density 
of nucleation centers is too large for a nucleus to be 
a rare event, and c) the growing probability is always 
negative, independent of the size of the nucleus. Nuclei 
which appeared immediately after the quench usually 
disappear again, almost none of them will grow to a 
domain in the polycrystal since the local order fluctuates 
strongly. Quantitative agreement with the heuristic 
Johnson-Mehl-Avrami-Kolmogorov description is ob- 
served and even for the deepest quench we did never 
observe an amorphous phase. We hope our work will 
stimulate theory and simulations to investigate such 
novel phenomena observed in phase transitions far from 
thermal equilibrium. 
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